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We derive distance priors from Planck first data release, and examine their impact on dark energy 
constraints from current observational data. We give the mean values and covariance matrix of 
{R,la,^bh ,ns}, which give an efficient summary of Planck data. The CMB shift parameters are 
R — \/Tlmh^r{zt), and la = 'Kr{z,)/rs{zt,), where z* is the redshift at the last scattering surface, 
and r{z*) and rs{zt) denote our comoving distance to z, and sound horizon at z* respectively. 

We find that Planck distance priors are significantly tighter than those from WMAP9. However, 
adding Planck distance priors does not lead to significantly improved dark energy constraints using 
current data, compared to adding WMAP9 distance priors. This is because Planck data appear 
to favor a higher matter density and lower Hubble constant, in tension with most of the other 
current cosmological data sets. Adding Planck distance priors to current data leads to a marginal 
inconsistency with a cosmological constant in a flat Universe. 
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I. INTRODUCTION 

Current observational data do not yet allow us to dif- 
ferentiate two likely explanations for the observed cosmic 
acceleration [l|, Q : dark energy, and the modification of 
general relativity. For recent reviews, see 04111- Cosmic 
acceleration is generally referred to as "dark energy" for 
convenience. 

There are three vigorously studied direct probes of 
dark energy. Type la supernovae (SNe la) probe the 
Hubble parameter H{z) (i.e., the expansion history of the 
Universe) via the measurement of luminosity distances 
to the SNe la [H, [1]. Galaxy clustering (GC) directly 
probes H{z) (and its integral form Da{z)) via the baryon 
acoustic oscillation (BAO) [12, llJ| measurements, and 
the growth rate fg{z) (i.e., the growth history of cosmic 
large scale structure) via redshift space distortion mea- 
surements. Weak lensing of galaxies probe a combination 
of the expansion history and growth history of the Uni- 
verse [11 El. 

While these direct probes of cosmic acceleration com- 
plement each other, each with its own set of system- 
atic uncertainties, they require the inclusion of cosmic 
microwave background (CMB) anisotropy data to help 
break the degeneracies among the dark energy and cos- 
mological parameters. This is because CMB data pro- 
vide the strongest constraints on cosmological parame- 
ters (see, e.g., [16]). 

Direct measurements of the Hubble constant (see, e.g., 
[17| ) also help break the degeneracy amongst the dark 
energy and cosmological parameters. Other data, e.g.. 
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gamma ray bursts |18l - l20l |. can help strengthen the dark 
energy constraints. 

In this paper, we derive distance priors from the Planck 
first data release, and examine their impact on dark en- 
ergy constraints from current observational data. 

We describe our method in Sec.H, present our results 
in Sec. HI, and conclude in Sec. IV. 



II. METHOD 

Our main goal is to derive distance priors from Planck 
data, and illustrate their impact on current observational 
data. For simplicity and clarity, we only use methods 
that give geometric constraints on dark energy in this 
paper. The constraints on the growth rate of cosmic 
large scale structure are degenerate with the geometric 
constraints (see, e.g., [2l|, [22]). We adopt a conservative 
approach by marginalizing over the growth constraints. 

Geometric constraints on dark energy are derived from 
the measurement of distances. The comoving distance to 
an object at redshift z is given by: 

riz) = cH^' \nk\-'/hmn[\nk\'/^r{z)], (1) 

n-A-[4^y Eiz)^Hiz)/Ho 

where sinn(x) — sin(x), x, sinh(x) for Qk < 0, ilj, = 0, 
and rife > respectively; and the expansion rate of the 
universe H(z) (i.e., the Hubble parameter) is given by 

H'{z)^(^l^' (2) 

= H^ [n„,{i + zf + a.(i + z)^ + f^fc(i + zf + nxXiz)] , 



where 17^ + fi^ + ^k + ^x = 1, and the dark energy 
density function X{z) is defined as 



X{z) 



^ Pxjz) 
~ PxiO)- 



(3) 



Note that Qr — f^m/(l + z^q) <^ ^m (with Zeq denoting 
the redshift at matter- rediation equahty), thus the fir 
term is usually omitted in dark energy studies at z ^ 
1000, since dark energy should only be important at late 
times. 



A. CMB data 

CMB data give us the comoving distance to the 
photon-decoupling surface r{z*), and the comoving sound 
horizon at photo-decoupling epoch rs{z:,) [23 ■ Wang & 
Mukherjee (2007) [2!] showed that the CMB shift pa- 
rameters 



Planck (post lenaing) + WMAP polarizatioi 
■ Gaussian with y,)=301-57. ff, =0-18 




FIG. 1: One-dimensional marginalized probability distributions 
of CMB shift parameter la derived from Planck, WMAP9, and 
WMAP7 data. 



R = ^VtrnHlr{z^)/c, 
la = 7rr(z,)/rs(z*), 



(4) 



together with oj^ = ftbh^, provide an efficient summary of 
CMB data as far as dark energy constraints go. This has 
been verified by [25| . Replacing wt, with z* gives identical 
constraints when the CMB distance priors are combined 
with other data [26[. Using ujt, instead of z*, is more 
appropriate in an MCMC analysis in which uji, is a base 
parameter. 

The comoving sound horizon at redshift z is given by 



rs{z) = 



.dt' 



= cH, 



dz' 



Eiz'Y 



cH- 



da' 



3{l + Rba')a''^E^{z') 



(5) 



where a is the cosmic scale factor, a = 1/(1 + z), 
and a^E'^(z) = il,n(a + Oeq) -I- ilfco^ -t- flxX{z)a'^, 
with acq = Sl,.ad/r2,„ = 1/(1 + Zcq), and z^^ = 
2.5 X 10^fl„ih'^{TcMB/2.7K)-^. The sound speed is 

Cs = l/y3{l + Rba), with R^a ^ 3pb/{4:pj), Rb = 
SlbOOnbh^iTcMB/^.TK)-"^. We take Tcmb = 2.7255. 

The redshift to the photon-decoupling surface, z.^,, is 
given by the fitting formula |27| : 



1048 [l + 0mi2A{nbh^) 



-0.738] 



[1 + giin^nh^'] 



where 



5i 

92 



0.0783 (l}fc/i^)-°-'38 

1 + 39.5 (fib/l2)0-763 

0.560 

l + 21.1(Ob/l2)l-81 



(6) 

(7) 
(8) 
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FIG. 2: One-dimensional marginalized probability distribution 
of CMB shift parameter R derived from Planck, WMAP9, and 
WMAP7 data. 



The redshift of the drag epoch Zd is well approximated 
by Hi 



Zd = 



1291 (^n/l) 



2\0.251 



1 + 0.659(a„/i2) 



0.828 



[l + bi{Qbh^) 



2\621 



(9) 



where 



6i = 0.313(r!™/i2)-"-4i9 [1 + 0.607(an/i')°-^^*] ,(10) 
62 = Q.238{n,nh^f-^^^. (11) 

FigsHlH] show the one-dimensional marginalized prob- 
ability distributions (pdf) of {la, R,i-^b,ns) from Planck 
[H, WMAP9 [3(|, and WMAP7 [16] data, for wx{z) = 
— 1 and one massive neutrino (with mass of 0.06eV), and 
without assuming a flat Universe. We have used the 



■ Planck (post lensing) + WMAP polarization 
-Gaussian with (nji") ^.0222Q. ct ,= .00030 
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FIG. 3: One-dimensional marginalized probability distribution of 
the dimensionless baryon density, a;;, = fli,h^, derived from Planck, 
WMAP9, and WMAP7 data. 
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FIG. 4: One-dimensional marginalized probability distributions of 
powerlaw index of primordial matter power spectrum, Us, derived 
from Planck, WMAP9, and WMAP7 data. 



Planck archiv data to obtain constraints on (la, R, uJb, ng) 
from Planck and WMAP9; this archiv does not include 
wx{z) 7^—1 cases without assuming a flat Universe. For- 
tunately, the constraints on (la, R,uJi,,ns) (including the 
pdf 's) are not sensitive to the assumption about dark en- 
ergy [3l|. Three sets of pdf's are shown in Figs[l]|4l 

(1) Planck-|-lensing-|-WP: Planck temperature data com- 
bined with Planck lensing, as well as WMAP polarization 
at low multipoles (/ < 23). This set represents the tight- 
est constraints from CMB data only at present. 

(2) WMAP: WMAP 9 year temperature and polarization 
data. 

(3) WMAP7: WMAP 7 year temperature and polariza- 
tion data. 

The Planck-flensing+WP pdf's in FigsHH well fitted 



by Gaussian distributions with the following means and 
standard deviations: 

(;„)= 301.57, a(y = 0.18 

(R) = 1.7407, CT(i?) = 0.0094 

(cjb) = 0.02228, cr(w6) = 0.00030 

(us) = 0.9662, cr(n,) = 0.0075. (12) 

The normalized covariance matrix of (/q, i?, Wbj"-s) is 



/ 1.0000 0.5250 

0.5250 1.0000 

-0.4235 -0.6925 

V -0.4475 -0.8240 



-0.4235 -0.4475 \ 

-0.6925 -0.8240 1 

1.0000 0.6109 

0.6109 1.0000 / 



(13) 



For comparison, we have also obtained the WMAP9 
constraints on {la,R,uJt,ns). The WMAP9 in FigsHH 
well fitted by Gaussian distributions with the following 
means and standard deviations: 

(la) == 302.02, cr(/a) ==0.66 

(i?) = 1.7327, cr(i?) = 0.0164 

(ujb) = 0.02260, CT(wfc) = 0.00053 

(n,) = 0.9719, a{n,) = 0.0143. (14) 

The normalized covariance matrix of {la, Rjt-^bTns) is 



1.0000 0.3883 -0.6089 -0.5391 

0.3883 1.0000 -0.5239 -0.6523 

-0.6089 -0.5239 1.0000 0.8563 

-0.5391 -0.6523 0.8563 1.0000 



(15) 



The WMAP7 constraints are from [31|. 

Since the primary GG data we use in this paper have 
been marginalized over Ug [32|, we should marginalize 
the GMB distance priors over Ug as well. This means 
dropping the 4th row and 4th column from the normal- 
ized covariance matrix of {la, R,uJb,ns), then obtain the 
covariance matrix for {la,R,uJb) as follows: 

C0YcMB{Pi,P]) = a-{Pi) Cr{pj) N0TT[iC0VcMB{Pi,Pj), 

(16) 
where i,j = 1, 2, 3. The rms variance (j{pi) and the nor- 
malized covariance matrix NormGovcMS are given by 
Eqs.([ni) and ([T3|) for Planck+lensing-fWP, and Egs. p^ 
and (dnj for WMAP9 respectively. 

CMB data are included in our analysis by adding the 
following term to the x^ of a given model with pi — 
la{z*), P2 = i?(2*),and p3 = Wfc: 



data 



XcMB = ^Pi [Covp^.^^ {pi , pj )] Apj , Api = Pi -p. 

(17) 
where pf"-*-"- are the mean from Eq. ([T^ and Eq. pi)) . 
and Go\^Qj^,jg is the inverse of the covariance matrix of 
[i?(z*),/a(-z*),Wf,] from Eq. (Uni) . Note that p4 = Us should 
be added if the constraints on Us are included in the GG 
data. 



B. Analysis of SN la Data 

SN la data give measurements of the luminosity dis- 
tance di^{z) through that of the distance modulus of each 

SN: 



fJ-o 



M = 5 log 



Mpc 



where CmBS,i, CmBCi, and Csc.i are the covariances be- 
tween TTiB, s, and C for the i-th SN. Note also that 
(T^ j includes a pecuHar velocity residual of 0.0005 (i.e., 
15bkm/s) added in quadrature [33| . 

The statistical and systematic covariance matrices, 
Cstat and Csys, are generally not diagonal [33], and are 



-|- 25 (18) given in the form: 



where m and M represent the apparent and absolute 
magnitude of a SN. The luminosity distance dL{z) = 
(1 -|- z)r{z), with the comoving distance r(z) given by 
Eq.©. 

Care must be taken in interpretating supernova dis- 
tances in an inhomogeneous Universe [38]. We use the 
compilation of SN la data by Conley et al. (2011) [33| . 
which include the SNe la from the first three years of the 
Supernova Legacy Survey (SNLS3), the largest homoge- 
neous SN la data set publicly available at present, and 
apply flux-averaging to reduce the systematic bias due to 
weak lensing maginification of SNe [34 36], as detailed in 
Wang, Chuang, & Mukherjee (2012) [31]- 

For a set of 472 SNe la, Conley et al. (2011) [33] 
give the apparent B magnitude, tub, and the covariance 
matrix for Am = vib — Timod, with 

m„iod = 51ogio Vl{z\s) ~ a{s - 1) + /3C + 7W, (19) 

where 'Dl{z\s) is the luminosity distance multiplied by 
Hq for a given set of cosmological parameters {s}, s is 
the stretch measure of the SN light curve shape, and C is 
the color measure for the SN. 7W is a nuisance parameter 
representing some combination of the absolute magni- 
tude of a fiducial SN la, M, and the Hubble constant 
Hq. Since the time dilation part of the observed lumi- 
nosity distance depends on the total redshift Zhei (special 
relativistic plus cosmological), we have |37J 



VL{z\s) = c-^Ho{l + z^,i)r{z\s), 



(20) 



where z and Zhoi are the CMB restframe and heliocentric 
redshifts of the SN. 

For a set of N SNc with correlated errors, we have 13311 



X 



Am^ 



Am 



(21) 



where Am is a vector with N components, and C is the 
N X N covariance matrix of the SNe la. 
Note that Attj is equivalent to A/iQ, since 

Am EE rns - mmod = [^b + a{s - 1) - /3C] - M. (22) 

The total covariance matrix is [33| 

C = Dgtat + Cstat + Csys, (23) 

with the diagonal part of the statistical uncertainty given 
byHl 



D 



stat,ii 



5(1 + zO 



lensing ' 
n 2 



host correction 



....2 I ^ 2 '2 I q2 '2 

_Zi(l + Zi/2)lnl0j '''' "'"■ ^'* 

+2aCmB.. - 2/3C„^c,» - 2a/3C,c,», (24) 



Cstat + Csys = Vo + a^Va+l3^Vb + 2aVoa-2mb~2a(3Vab. 

(25) 
where Vq, K, Vb, Voa, Vob, and Kb are ma- 
trices given by the SNLS data archive at 
https://tspace.library.utoronto.ca/handle/1807/24512/. 
Cstat includes the uncertainty in the SN model. Cgys 
includes the uncertainty in the zero point. Note that 
Cstat and Csys do not depend on Ai, since the relative 
distance moduli are independent of the value of Ai [33] . 

We refer the reader to Conley et al. (2011) [33] for 
a detailed discussion of the origins of the statistical and 
systematic errors. As an example, we note that the cor- 
relation of errors on different SNe arises from a statistical 
uncertainty in the zero point of one passband, e.g., tm- 
This directly affects all SNe with tm measurements due 
to K-corrections (restframe B to vm), and indirectly af- 
fects even the SNe without tm measurements through 
the empirical SN models by changing the templates and 
the measured color-luminosity relationship. 

For x^ statistics using MCMC or a grid of parameters, 
here are the steps in flux-averaging [3l|: 

(1) Convert the distance modulus of SNe la into 
"fluxes" , 



F{zi) = 10-(A'o""'(^i)-25)/2.5 



^data, 



Zl) 



Mpc 



(26) 



(2) For a given set of cosmological parameters {s}, ob- 
tain "absolute luminosities", {£(z/)}, by removing the 
redshift dependence of the "fluxes", i.e., 

C{zi) = dl{zi\s)F{zi). (27) 

(3) Flux-average the "absolute luminosities" {£;} in 
each redshift bin i to obtain -^ 'C > : 



Ni 



£' 



N, 



^E^K^F^ 



1=1 






Ai) 



(28) 



1=1 



(4) Place C at the mean redshift Zi of the i-th redshift 
bin, now the binned flux is 

Fi-z,)=T/dli-z,\s). (29) 

(5) Compute the covariance matrix of F{zi) and F{zj): 

Cov[F{zi),Fizj)] (30) 

In 10/2.5 



N,Nj 



dL{zi\s)dL{zj\s) 



^ ^ /:(z«)/:(z(i))(A/.r^(z«)AMr^(zW)) 



1=1 m=l 



where (A^q'^''*(z}'^)A/^o'*'''(z2')) is the covariance of the 
measured distance moduh of the Z-th SN la in the «-th 
redshift bin, and the m-th SN la in the j'-th redshift bin. 
£(z) is defined by Eqs-dlSl) and (E^. 

(6) For the flux-averaged data, {F(zi)}, compute 

X'^Y. ^^(^^) Gov-' [Fi-z,),Fiz,)] AF(z,) (3f ) 



where 



AFCz,) = F{-z,)-FP{-z,\s), 



(32) 



with FP(z,|s) = {dL{z\s)/Mpcy\ 

For the sample of SNe we use in this study, 
we flux-averaged the SNe with dz ~ 0.07, to 
ensure that all redshift bins contain at least one 
SN. Our SN flux-averaging code is available at 
http://www.nhn.ou.edu/~wang/SNcode/ 



C. Galaxy Clustering Data 

For GC data, we use the measurements oi H{z)rs {zd)/c 
and DA{z)/rs{zd) (where H{z) is the Hubble parame- 
ter, Da{z) is the angular diameter distance, and rg^zj) 
is the sound horizon at the drag epoch) from the two- 
dimensional two-point correlation function measured at 
z=0.35 [34I and z=0.57 [33|- The z = 0.35 measurement 
was made by Chuang & Wang (2012) [32| using a sample 
of the SDSS DR7 Luminous Red Galaxies (LRGs). The 
z = 0.57 measurement was made by Chuang et al. (2013) 
[H using the CMASS galaxy sample from BOSS. 

Using the two-dimensional two-point correlation func- 
tion of SDSS DR7 in the scale range of 40-120 Mpc//i, 
Chuang & Wang (2012) [H found that 

H{z^O.Zb)rs{zd)/c = 0.0434 ±0.0018 
DA{z^Q.ib)/rs{zd) = 6.60 ±0.26 

r = 0.0604 (33) 

where r is the normalized correlation coefficient between 
H{z = 0.35)rs(Z(i)/c and Da{z = Q.i?>)/rs{zd)j and 
rs{zd) is the sound horizon at the drag epoch (given by 
Eqs.® and ^. 

In a similar analysis using the CMASS galaxy sample 
from BOSS, Chuang et al. (2013) found that 

H{z = Q.^l)rs{zd)/c = 0.0454 ± 0.0031 
DA{z = Q.^l)lrs{zd) = 8.95 ±0.27 

r = 0.4874 (34) 

We marginalize over the growth rate measurement made 
by Chuang et al. 2013 [39] for a conservative approach. 
GC data are included in our analysis by adding ^^(j = 
Xgci +XGC2' with ZGci = 0.35 and zgc2 = 0.57, to the 
X^ of a given model. Note that 

X^c« = Ap^Cgc ^V^ , Pj )] Apj , Ap, = p, - vf , 

(35) 
where pi = H{zGCi)rsizd)/c and P2 = DA{zGCi/rs{zd), 
with i — 1,2. 



D. Gammay-ray Burst Data 

We add gammay-ray burst (GRB) data to our analy- 
sis, since these are complementary in redshift range to the 
SN la data. We use GRB data in the form of the model- 
independent GRB distance measurements from Wang 
(2008c) [40|, which were derived from the data of 69 
GRBs with 0.17 < z < 6.6 from Schaefer (2007) [4l|i. 

The GRB distance measurements are given in terms of 

m 



rp{zi) ^ 



(^) 



rp(0.17)^ 



rpiz) 



_(1 



-zy/' Ho , , 
z en 



(36) 
where r{z) is the comoving distance at z. 

The GRB data are included in our analysis by adding 
the following term to the x^ of a given model: 



Xgrb = [Arp(zi)] • (Covq^s) . • [Arp{zj)] 



Ar^(z,) = 7v^^'^(z,)-7>(z,), 



(37) 



where rp{z) is defined by Eq. (j36p . The covariance matrix 
is given by 

iCovGB.B)ij = cr(r>(2;0)cr(r>(zj)) (CovGfls),^- , (38) 

where Goygrb is the normalized covariance matrix from 
Table 3 of Wang (2008c) "ic^, and 



f — / ^^ + 



— — /^\ ^ tt;— /" .^^da'ta. 



CT{rp{zi)) = (7 {rp{zi)) , a rp{z) > rp{z) 
<j{r^{z,)) = a{r^{z,))-, ii r^{z) < r^{zf-'M39) 

where a {r^{zi)) and a (¥^{zi))^ are the 68% C.L. errors 
given in Table 2 of Wang (2008c) |40|. 



E. Dark energy parametrization 

Since we are ignorant of the true nature of dark energy, 
it is useful to measure the dark energy density function 
X{z) = pxiz)/px{0) as a free function of redshift [46l - 
l48l | . This has the advantage of allowing dark energy mod- 
els in which px{z) becomes negative in the future, e.g., 
the "Big Crunch" models [4^, [5^ , which are precluded 
if we parametrize dark energy with an equation of state 
wx{z) \0i. 

Here we parametrize X{z) by cubic-splining its values 
at z = 1/3, 2/3, and 1.0, and assume that X{z > 1) = 
X{z = 1). For simplicity of notation, we define X0.33 = 
X{z = 1/3), Xo.67 = X{z = 2/3), and Xi.o = X{z = 
1). Fixing X{z > 1) reflects the limit of current data. 



^ The proper calibration of GRBs is an active area of research. 
For recent studies on the impact of detector thresholds, spectral 
analysis, and unknown selection effects, see, e.g., [42l445|l . 



and avoids making assumptions about early dark energy 
that can be propagated into artificial constraints on dark 
energy at low z [2J, |43| . 

For comparison with the work of others, we also con- 
sider a dark energy equation of state linear in the cosmic 
scale factor a, wx(q) = wq ^ (\ — a)wa ^51]. A related 
parametrization is [2 a ) 



SNLS3{nux-averaged)+HST(Rn)+GRB+BA0(CW2)+: Planck+lensing+WP (solid); WMAP (dotted) 



wx{.z) = wo(3a - 2) + 3wo.5(l - a), 



(40) 



where wq.s = wx{z = 0.5). Wang (2008b) [26] showed 
that (iuq, iL'o.s) are much less correlated than (wq, Wq), 
thus are a better set of parameters to use. We find 
that (wq, Wo. 5) converge much faster than (wq^ Wa) in 
a Markov Chain Monte Carlo (MCMC) likelihood anal- 
ysis for the same data. 




III. RESULTS 



We perform a MCMC likelihood analysis [52| to obtain 
0(10®) samples for each set of results presented in this 
paper. We assume flat priors for all the parameters, and 
allow ranges of the parameters wide enough such that fur- 
ther increasing the allowed ranges has no impact on the 
results. The chains typically have worst e- values, defined 
to be the variance(mean)/mean(variance) of 1/2 chains, 
much smaller than 0.005, indicating convergence. The 
chains are subsequently appropriately thinned to ensure 
independent samples. 

In addition to the SN la, CMB, GC, and GRB data 
discussed in SecHIl we impose a prior of Hq = 73.8 ± 
2.4kms~^Mpc~^, from the HST measurements by Riess 
et al. (2011) [IZ]- 

We do not assume a flat universe. In addition to the 
dark energy parameters described in Sec lIIE| we also 
constrain cosmological parameters (r2,„, fi^, h, ujb), where 
uji, = Vli,h? . In addition, we marginalize over the SN 
la nuisance parameters {a,/3,A^}. We only use flux- 
averaged SN la data (with dz = 0.07), as flux-averaging 
reduces the impact of systematic uncertainties on dark 
energy and cosmological parameter constraints [3l| . 

We will present results for dark energy density at z = 
1/3, 2/3, and 1, as well as (wo,Wa) and (wq, W0.5): and a 
constant dark energy equation of state w. 



FIG. 5: The marginalized probability distribu- 

tions for {XQ,2,-i,XQS7,X\.o,Q.m,^k,h,u)i,,a,P,M}, for 

SNe-|-Ho-|-GRB-|-GC(CW12), combined with Planck-|-lensing-(-WP 
(solid) and WMAP9 (dotted) data respectively. 
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FIG. 6: The joint 68% and 95% confidence contours for param- 
eters of interest for SNe-|-_ffo-f GRB+GC(CW12), combined with 
Planck+lensing-(-WP (solid) and WMAP9 (dotted) data respec- 
tively. The thin solid and dotted contours also include GC data 
from BOSS. 



Constraints on dark energy density function 

X{z) 



A. 



Figsl5][5] summarize our constraints on A(z) 
parametrized by its value aX z — 1/3, 2/3, and 1. 
Planck data give very similar results as WMAP9 data 
on X{z), even although Planck data favor higher Vim- 
However, note that adding Planck priors leads to a 
marginal inconsistency with a cosmological constant in 
a flat Universe (see bottom right panel in FiglB]). 



Adding BOSS data has a more significant impact: it 
shifts the value of Ag.e? away from 1 at ~ 1.5(7 (see Figl6] 
and FigjS]) , independent of cosmic curvature. 



B. Constraints on a linear dark energy equation of 

state 

We have studied the constraints on both {wo,Wa) and 
(u'o,wo.5), as these have different base parameters (as- 
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FIG. 7: The 68% and 95% confidence constraints 

for X{z) (cubic-splined from {^o.33i ^0.67i^i.o})i for 
SNe+Ho+GRB+GC(CW12), combined with Planck+lensing+WP 
(soUd) and WMAP9 (dotted) data respectively. 



SNe(f.a.)+GC(CW1 2)+Hj+GRB+ 



SNe(f.a.)+GC(CW12)+H +GRB 
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FIG. 9: The joint 68% and 95% confidence contours 
for {wo,Wa) (left panel) and {wo,wo.5) (right panel), for 
SNe+^0+GRB+GC(CW12), combined with Planck+lensing+WP 
(solid) and WMAP9 (dotted) data respectively. 
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FIG. 8: The 68% and 95% confidence constraints 

for X{z) (cubic-splined from {^0.33, -''^o.eTj^i.o}), for 
SNe+Ho+GRB+GC(CW12+BOSS), combined with 

Planck+lensing+WP (solid) and WMAP9 (dotted) data re- 
spectively. 



FIG. 10: The joint 68% and 95% confidence contours 
for (Oj;,uia) (left panel) and (f2fc,uio.5) (right panel), for 
SNe+Ho+GRB+GC(CW12), combined with Planck+lensing+WP 
(solid) and WMAP9 (dotted) data respectively. 



sumed to have flat priors). In order to compare with 
previous work, and to display the impact of replacing 
WMAP9 priors with Planck+lensing+WP priors, we do 
not include GC data from BOSS in this comparison. 

Fig[S] shows the joint 68% and 95% confidence con- 
tours for {wQ,Wa) (left panel) and (^0,^0.5) (right 
panel), for SNe-HiJo+GRB+GC(CW12), combined with 
Planck+lensing+WP (sohd) and WMAP9 (dotted) data 
respectively. FigHU] shows the corresponding joint 68% 
and 95% confidence contours for {Qk,Wa) (left panel) and 
(i7fe,wo.5) (right panel), for the same data and with the 
same line types. Fig jS] indicates that Planck priors do not 
have a significant impact on the constraints on a linear 
dark energy equation state; this means that the DETF 
FoM remains approximately the same compared to that 
found by [3l| using WMAP7 priors (using WMAP9 priors 
gives similar results as using WMAP7 priors). It is inter- 
esting to note that the right panel of FigHO] shows that 
the combined data with Planck priors rule out W0.5 = — 1 
and a flat Universe at ^ 1.5a. Planck data favor a small 
but positive 51^ (i.e., a slightly open Universe). 
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FIG. 11: The joint 68% and 95% confidence contours for (Qj;,to) 
(left panel) and (Qm, w) (right panel), for SNe+_ffoi combined with 
Planck+lensing+WP (solid) and WMAP9 (dotted) data respec- 
tively. 



Constraints on a constant dark energy equation 
of state 



In order to understand better the difference between 
the Planck+lensing+WP priors and the WMAP9 priors, 
we now study a constant dark energy equation of state, 
for the minimal combination of SNe+i/o data with the 
CMB priors. 

Fig[TT] shows the joint 68% and 95% confidence con- 
tours for {rik , w) (left panel) and {flm , w) (right panel) , 
for SNe+iJo, combined with Planck+lensing+WP (solid) 
and WMAP9 (dotted) data respectively. Here, we see 
even more clearly the trend of Planck+lensing+WP pri- 
ors favoring a small but postive flk (a slightly open Uni- 
verse), and a somewhat higher fi™, compared to WMAP9 
priors. Again, we find that adding Planck priors to other 
data leads to a marginal inconsistency with a cosmolog- 
ical constant in a flat Universe. 



IV. DISCUSSION AND SUMMARY 

We have derived the distance priors from Planck first 
data release, in the form of the mean values and covari- 
ance matrix of {i?, /aj ^b^^,'^s}, which give an efficient 
summary of Planck data in the context of dark energy 
constraints. We have used constraints on {i?, Z^, Jlfc/i^} 
from Planck data in combination with other data to 
probe dark energy in a conservative geometric approach. 

We have considered three different dark energy 
parametrizations: (1) Dark energy density X{z) — 
Px{z)/ px(0), parametrized by splining its values at z = 
1/3, 2/3, and 1.0. (2) Dark energy equation of state lin- 
ear in the cosmic scale factor, wxiz) = wq + Wa(l — a), 
as well as the alternative parametrization with much less 
correlated parameters uiQ andwo.5, wxiz) — wo(3a — 2) + 
3^0.5(1 — a)- (3) Constant dark energy equation of state 
wx{z) = w. 

In addition to C MB p riors, we used SNe compiled by 
Conley et al. (2011) (Sj), flux-averaged to reduce system- 



atic errors, the Hubble constant measurement by Riess 
et al. (2011) [13], the GRB data as summarized in Wang 
(2008) m, and GC data from SDSS DR7 derived by 
Chuang & Wang (2012) [32|. We have chosen to use 
the same data as in [3l|, except for the CMB priors 
(Planck+lensing+WP and WMAP9, versus WMAP7), 
in order to compare the impact of the various CMB pri- 
ors. For completeness, we have added the GC data from 
BOSS derived by Chuang et al. (2013) [33| in constrain- 
ing X{z). 

We find that when dark energy density is allowed to 
be a free function, current data (excluding GC data from 
BOSS, with either Planck+lensing+WP or WMAP pri- 
ors) are fully consistent with a cosmological constant and 
a flat Universe at 95%, but deviate from a cosmological 
constant in a flat Universe at ~ 68% (see Figl?]). The 
addition of BOSS data leads to a deviation from a cos- 
mological constant at the redshift near that of the BOSS 
data at > 68% confidence independent of cosmic curva- 
ture (see Fig[5]). In general, adding Planck+lensing+WP 
priors leads to a preference for a small positive flk (i-e., 
a slightly open Universe) for a cosmological constant, or 
a flat Universe with dark energy deviating from a cos- 
mological constant, compared to adding WMAP9 priors 
(see FigH). 

When dark energy equation of state is assumed to 
be a linear function in the cosmic scale factor a(i), the 
dark energy constraints depend on the base parameters 
used, the highly correlated {wo^Wa}, or the much less 
correlated {^0,^0. 5}- The constraints on {wo^Wa} are 
consistent with a cosmological constant and a flat Uni- 
verse for both Planck+lensing+WP and WMAP9 priors, 
while that on {wq, wo.5} are marginally inconsistent with 
a cosmological constant in a flat Universe (see FiglHland 
FigfTO)) for Planck+lensing+WP priors, similar to our 
findings in the X(z) case. 

We find that the above trend in the 
Planck+lensing+WP versus WMAP9 comparison 
becomes more pronounced when we assume a con- 
stant dark energy equation of state. Here we only 
combine CMB priors with a minimal set of other 
data, SNe and Hq, to highlight the difference between 
Planck+lcnsing+WP and WMAP9 priors. We find 
that adding the Planck+lcnsing+WP priors to the SNe 
and Hq measurement leads to an inconsistency with a 
cosmological constant in a flat Universe (see Fig fTTj) at 
> 68% confidence. 

To conclude, we find that Planck distance priors are 
significantly tighter than those from WMAP9 (see Figs[l]- 
[3]). However, adding Planck distance priors does not lead 
to significantly improved dark energy constraints using 
current data, compared to adding WMAP9 distance pri- 
ors. This is because Planck data appear to favor a higher 
matter density and lower Hubble constant [29| , in tension 
with most of the other current cosmological data sets. 

In order to understand the nature of dark energy, we 
will need to improve our understanding of the systematic 
uncertainties of all data used. Future dark energy mea- 



surements from space |53l - l56| that minimize systematic 
uncertainties by design will enable us to make dramatic 
progress in our quest to shed light on dark energy. 
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